Stochastic versus dynamic approach to Levy statistics in the presence of an external 

perturbation 
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We study the influence of a dissipation process on diffu- 
sion dynamics triggered by slow fluctuations. We study both 
strong- and weak-friction regime. When the latter regime ap- 
plies, the system is attracted by the basin of either Gauss 
or Levy statistics according to whether the fluctuation cor- 
relation function is integrable or not. We analyze with a 
numerical calculation the border between the two basins of 
attraction. 



I. INTRODUCTION 

There is a growing interest in literature |l],|^] on the 
physical manifestation of Levy diffusion. This interesting 
subject can be dealt with in two distinct ways. The first 
rests on the assumption that there exist in nature Levy 
fluctuations, namely stochastic processes of the Levy sta- 
ble form Hf|. We refer to this approach as stochastic. 
There are many examples of this attitude, and we limit 
ourselves to quoting a sample of the early work based on 
this view . One of the most recent examples of this 
kind of approach is given by the work of Ref. |8[ . 

The second way of dealing with Levy fluctuations is 
based on the assumption that these processes have a 
Hamiltonian foundation. We refer to this kind of ap- 
proach as dynamic. An outstanding champion of this 
view is Zaslavsky and the interested reader is referred 
to his recent book || for a transparent illustration of 
this perspective. According to Zaslavsky's analysis, the 
Hamiltonian dynamics yielding Levy diffusion is charac- 
terized by specific weak-chaos properties. For calcula- 
tion purposes these properties can be mimicked by either 
maps [fU)|-[i"2"[ of the same kind as those used to study in- 
termittency jl3| ] or a suitable non-linear transformation 
of a random noise generator [|14|. Therefore we shall refer 
ourselves to all these papers [|To| [l^ , |T^ as examples of a 
dynamic approach to Levy diffusion not to speak of those 
papers explicitly using Hamiltonian systems as diffusion 
generators |l5| , |l6| ]. 

As far as the case of free diffusion is concerned, the 
dynamic and stochastic approach yield almost indistin- 
guishable results. The main purpose of this letter is to 
show that when the interesting case of a perturbation 



is considered, the dynamic approach and the stochastic 
approach can lead to totally different predictions. A rec- 
onciliation of the two approaches can only be obtained in 
the limiting case of external perturbations of extremely 
weak intensity. 



II. STOCHASTIC APPROACH 

We consider a free motion of a point particle subject to 
a Markovian stochastic noise 77(f). This means that for 
any infinitesimal time interval dt there is a momentum 
change per unit of mass given by: 



dv(t) = dr](t) 



(1) 



where v(t) is the velocity of the particle. In the treatment 
of this letter the forcing term will be either a Wiener or 
a Levy process. This means that the probability density 
for the fluctuation 77(f) to make a jump by the intensity 77 
in the time interval f , P(r), t) , defined through the inverse 
Fourier transform of its characteristic function, reads: 
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in the case of a Wiener process, and: 
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in the case of a symmetric Levy process. The Levy pro- 
cess is characterized by a positive parameter b, determin- 
ing the width of the distribution and by the coefficient a 
determining the specific Levy statistics of the system. 

As discussed in Ref. [H, the dynamic approach to 
Levy processes has to be restricted to the range: 1 < 
a < 2, while, as well known J|Jl7]], the full range of Levy 
processes is: < a < 2. Thus, to make a comparison 
between the stochastic and the dynamic approach we are 
forced to restrain our analysis to the restricted interval 
of a values: 1 < a < 2. In this interval we do not 
have available any known analytical expression for the 
density distribution P(r),t). However, we know [jlTj that 
the distribution is bell-shaped and that in the asymp- 
totic limit of large 77 is characterized by inverse power 
law tails, proportional to f/|77|( 1+a ). As far as the width 
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of the distribution is concerned, it cannot be defined by 
the variance, which is infinite. In a loose sense the width 
of the distribution is determined by the parameter b. It is 
evident that Eq.(|l|) makes the distribution P(v,t) iden- 
tical to P(r), t) 

Let us consider the perturbation of the free diffusion 
generated by Eq. (jlj) by means of: 



dv(t) = K(v) dt + dr/(t) 



(4) 



where K(v) is a generic perturbation. For simplicity 
we restrict our analysis to the case of linear damping, 
thereby replacing K{v) in Eq.(|J) with — Xv(t): 



dv(t) = -Xv(t)dt + dT)(t) 



(5) 



When the fluctuation dr]{t) is the Wiener process of Eq. 
(g), Eq.(||) becomes the ordinary Ornstein-Uhlenbeck 
process, which yields the Fokker-Planck equation and 
with it the equilibrium Gaussian density, whose Fourier 



transform is e 2 k ; with variance: 
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If r](t) is the Levy process of Eq.(||), the Fourier transform 
of the equilibrium distribution is e~ hA l fe l°, with the same 
exponent a as that of free diffusion. The parameter b\ 
(see |Mp|) is defined by: 
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(7) 



We note that b\ and er? share the same structure. In both 
cases the effect of the friction term is that of quench- 
ing the free diffusion so as to generate a time inde- 
pendent, or equilibrium, distribution. This means that 
Iim t -»ooi , (»,*) = V(v). 



III. DYNAMIC APPROACH 

Both the Wiener and Levy noise illustrated in Section 
II are mathematical abstractions with a limit of validity. 
The discrepancy between this mathematical abstraction 
and physical reality can become significant as a result 
of an external perturbation forcing the dynamics of the 
system to produce physical effects stemming from the 
time scale where the mathematical idealization departs 
from physical reality. Here we focus our attention on the 
case where the physical reality is assumed to be properly 
described by the dichotomous fluctuation £ used in earlier 
work for a dynamic derivation of Levy processes [fl4[ . 

We assign to this dichotomous variable the values W 
and —W. According to Ref. jL4j the statistics of the dif- 
fusion process generated by this variable are determined 
by both its dichotomous nature and its correlation func- 
tion <&{(£)• We assign to this correlation function the 
form: 
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(8) 



According to the theoretical analysis of Ref. 14 jig], the 
dichotomous nature of this fluctuation makes especially 
relevant the physical meaning of the function ip^(t) de- 
fined by: 



{pr + t)0H 



This function is the distribution of sojourn times in one of 
the two equiprobable states of the dichotomous variable 
£. The parameter T denotes the mean waiting sojourn 
time. 

It is convenient to illustrate some aspects of the free 
diffusion process generated by this fluctuation, namely, 
the case described by: 

dv(t) = f(t) dt . (10) 

A realization of this process of free diffusion is: 
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«(t)=M*-*n-i) + Z>n 
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(11) 



where t^ < t < t ( °\ + t„ 



n-l 
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Here the r' k s denote the time durations of the sojourns in 
the accelerating states occurring prior to time t, and the 
£' k s denote cither the value W or the value — W of the 
dichotomous fluctuation. The distribution of the sojourn 
times is given by Eq. (^). Thus we easily derive the 
distribution of the velocity jumps Avk — CfcTfc, which 
turns out to be: 
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(12) 



This distribution is of central importance to predict the 
statistics of the resulting diffusion process. If (3 > 1 the 
second moment of this distribution is finite. Thus we are 
allowed to apply the conventional central limit theorem, 
thereby recovering a diffusion process, indistinguishable 
in the long-time limit from those triggered by the Wiener 
fluctuations, with variance per unit of time given by: 



2f3W 2 T 
/3-1 



(13) 



We say J!?]] that the diffusion process is attracted by 
the basin of Gauss statistics. If < 8 < 1 the sec- 
ond moment of the distribution of Eq.(|l2J) is not finite, 
and we are therefore forced to use the Levy-Khintchine- 
Gnedenko theorem |l7j (that is, the extension of the cen- 
tral limit theorem for random variables with no finite 
variance). According to Ref. [17] and to the result of 
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the Appendix, this condition yields a Levy process with 
parameters b and a given by: 



b = W 



(/3T^sin(|/3)r(l-/3) (14) 



and 



a = (3 + 1 



(15) 



In this case we say that the diffusion process is attracted 
by the basin of Levy statistics. 

We see that in the long-time limit the diffusion gener- 
ated by £ is essentially indistinguishable from that gener- 
ated by the Levy fluctuation r\. What about the influence 
of an external perturbation on these two distinct fluctu- 
ations? One might be tempted to make the conjecture 
that both processes result in the same equilibrium dis- 
tribution. However, for this conjecture to be correct it 
would be necessary to apply first the generalized central 
limit theorem and, afterwards, the external perturbation. 
In general, the result is different from what one would 
obtain reversing this order, namely, applying the pertur- 
bation first, and the time asymptotic analysis next. Thus 
the naive prediction that both fluctuations produce the 
same equilibrium distribution, in general is wrong. 

To substantiate our view, let us consider the dynamic 
counterpart of Eq. (||) : 



- V (t) + \ V (t) = m 



(16) 



It is easy to solve Eq.(16). In a single motion event we 
have: 



/ x ±W , 

v(t) = —(!-< 
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(17) 



with 



t {0) <t<t {0) 
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and 



fc=0 



Here v(t) is forced to remain in the interval 
[— W/X, +W/X}. The computer simulation of this letter 
is based, as done in Ref. |L9|, on generating, by means 
of a suitable non-linear deformation of the random num- 
bers with a uniform distribution in the interval [0,1], a 
sequence of sojourn times with the density of Eq. , 
and with another random number generator, equivalent 
to tossing a coin, a random sequence of velocity signs. 
For each trajectory realization we set the initial con- 
dition v(0) = and the trajectory is observed till to 
a fixed stop time, and subsequently recorded in a bin. 
Each trial was repeated 10 4 times. In all the numerical 
calculations behind the results illustrated in the enclosed 
figures, we assume that equilibrium is reached after a 
time > 20/ A. The following subsections are devoted to 
illustrating the results obtained in three distinct physical 



conditions. Note that to illustrate the equilibrium dis- 
tribution shape corresponding to different values of A we 
use a variable z, which is the rescaled velocity obtained 
by multiplying the original velocity v by the factor W/X. 
Thus the rescaled velocity ranges from z = — 1 to 2 = 1. 



A. Large perturbation 

In the strong damping case, say when XT (3 > 1, for 
any of the two accelerating states equilibrium is reached 
before the state comes to its end. This means that the 
preferred velocity values will be either W/X or —W/X. 
This also means that the equilibrium distribution is ex- 
pected to be U-shaped, as confirmed by the curve of Fig. 
|l| corresponding to A = 0.05. This damping-induced U- 
shaped distribution is generated by system's dynamics 
regardless of whether the corresponding free diffusion is 
located in the Levy or the Gauss basin of attraction. 



B. Weak perturbation 

We have seen that in the case of very large friction the 
distribution has a distinct U-shaped form: This means 
two peaks enclosing an almost empty central region. 
From Fig. g we see also that weakening the perturba- 
tion intensity has first the effect of populating this empty 
central region with a uniform distribution (see the curve 
corresponding to A = 5 • 10~ 3 ). A further decrease of the 
perturbation intensity (A < 4 • 10 -4 ) makes a bell-shaped 
Levy distribution of decreasing width and increasing am- 
plitude appear in this central part of the equilibrium dis- 
tribution. 

Fig. U is devoted to illustrating the effect that chang- 
ing (3 has on the distribution shape, in a case of weak 
friction, namely, with A = 10~ 5 . Note that in the Gauss 
basin of attraction, with (3 > 1, the distributions are 
Gaussian functions with no algebraic tails, whereas in 
the Levy basin of attraction, with (3 < 1, the equilibrium 
distribution is markedly characterized by the emergence 
of slow tails. The transition from the one to the other 
condition takes place at (3 = 1. 



C. Transition from the Levy to the Gauss basin of 
attraction 

As pointed out in subsection 3 B, in the weak- friction 
regime, the statistics of the equilibrium distribution are 
determined by the values of the parameter (3. According 
to whether system's dynamics belongs to the Gauss or 
the the Levy basin of attraction, we are led to adopt dif- 
ferent criteria to determine the distribution width using 
the stochastic analysis, i.e. the variance er? or the pa- 
rameter b\. In the former case, (3 > 1, the variance a\ is 
obtained using Eq.(0) and Eq.fO). This yields: 
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FIG. 1. Equilibrium probability densities for different values of the friction A at (3 — 0.6. The Levy shape emerges upon 
friction decrease. 




FIG. 2. Equilibrium probability densities for different values of (3 at A = 10 5 . The distribution change shape from Levy to 
Gauss one when (3 crosses the critical value (3 = 1. 
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FIG. 3. The distribution widths as a function of (3. The curves are parametrized by the friction A. For (3 < 1 the ordinates 
refer to b\. For (3 > 1 the ordinates refer to a\. The lines illustrate the theoretical prediction according to the stochastic 
analysis. The points are the numerical results. 
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FIG. 4. The population of the region 1 > \z\ > 0.8 as a function of /3, for different values of the friction A. In the small 
friction condition the population quickly drops to zero with (3 moving from f3 < 1 to (3 > 1. 
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(18) 



thereby resulting in a divergence for f3 — > 1 + . On the 
r.h.s. of Fig. || we establish a comparison between this 
theoretical prediction and the corresponding numerical 
result. For (3 >> 1 we find an excellent agreement be- 
tween theory and numerical simulation. We note however 
that the numerical calculations do not produce any diver- 
gence at the transition from the Gauss to the Levy basin 
of attraction. In the latter case, when j3 < 1, the ex- 
pected Levy equilibrium distribution (see |Q) should be 
characterized by the value of the parameter b\ obtained 
using Eqs.(g), @, (|l|). This value is: 



W{0TWY 

W+T) 



in T(l - 0) , 



(19) 



which results again in a divergence at the transition from 
the Levy to the Gauss basin of attraction. In the left 
side of Fig.|| the comparison between theoretical predic- 
tion and numerical calculation shows a good agreement 
around (3 ~ 0.5. 

Although the numerical calculation does not result in 
any divergence, we see that decreasing the friction inten- 
sity has the effect of improving the agreement between 
theoretical prediction and numerical finding. The mea- 
sured values of the parameters 6a and <j\ becomes in- 
creasingly larger at the transition from one basin of at- 
traction to the other, and the overall behavior becomes 
increasingly similar to that of a phase transition. 

Fig. |J is devoted to illustrating the peak contributions 
to the equilibrium distribution. We have evaluated nu- 
merically the amount of population for 0.8 < \z\ < 1. We 
see that decreasing A has the effect depleting this region 
if (3 > 1, while a significant amount of population is left 
in this region if (3 < 1. This numerical result shows that 
at extremely small friction values, moving from (3 < 1 
to j3 > 1 has the significant effect of making the peak 
intensity drop to zero, even if the phase-transition char- 
acter of the passage from [3 < 1 to (3 > 1 is characterized 
by rare intense fluctuations and, consequently, by a large 
numerical error. 



IV. CONCLUDING REMARKS 

As its main contribution, this paper sheds light into 
the difference between the conventional stochastic ap- 
proach and the dynamic approach. Of some relevance 
is the numerical method adopted. Rather than using as 
stochastic generators intermittent maps we have founded 
our numerical treatment on the random generation of the 
waiting time distribution, a fact that allowed us to set- 
tle numerically problems that would have implied other- 
wise hard numerical difficulties. Setting apart the case 
of strong friction where the stochastic and the dynamic 
method yield strikingly different results, the numerical 



method used made it possible to settle two much more 
delicate problems: 

(i) We have pointed out the residual differences be- 
tween the two method surviving in the extremely weak 
friction regime. 

(ii) We have made accessible to numerical investigation 
the delicate issue of the transition occurring at j3 = 1 
from the Gaussian ((3 > 1) to the Levy (/? < 1) statistics. 

It has to be pointed out that this transition region is 
not yet well understood Q and further research work 
along the lines of this letter might serve the interesting 
purpose of establishing whether in the limiting case of 
very small friction the equilibrium distribution yields the 
divergencies predicted by Eq. ( |l8| ) and Eq.(19) or a finite 
values, as suggested by the recent theoretical analysis of 
Ref. pc[| . We see from Fig. ||that the numerical analysis 
yields finite rather than divergent widths at [3 = 1, but 
at the moment it is not yet clear if this is due to T > 0, 
as argued in Ref. [20], or to the adoption of not yet 
sufficiently small values of A. 



APPENDIX: 

In this Appendix we show how to determine the pa- 
rameter b defining the width of the diffusion process (||) 
generated, in the long-time limit, by the dichotomous 
variable of Eq.([l0|). To realize this goal we use some fun- 
damental theorems established by Levy, Kintchine and 
Gnedenko, whose detailed demonstration can be found 
in Ref. Q. 

The central issue is the assessment of the limit dis- 
tribution of the normalized sum of the independent and 
identically distributed random variables £i, £2, • • • , C« : 



Cl + C2 + ■ ■ ■ + Cn 



(Al) 



where A n and B n are suitable normalization constants. 
To establish this limit condition we rest on the following 
three properties. 

First, we define what a stable distribution is all about. 
A distribution V{x) is stable if given the arbitrary real 
numbers a\ > 0, b%, a 2 > 0, b 2 there exist the numbers 
03 > and &3 such that the equality 

V{aix + 61) * V(a 2 x + b 2 ) = V{a 3 x + 63) (A2) 

holds. Here * denotes the convolution product. 

Second, we use a theorem by Kintchine and Levy |T^ ] 
which establishes that the stable distributions are t hose, 
and only those, corresponding to the sum of Eq. ( |Alj ) con- 
verging to a finite limit as n — > 00. 

Third, using another theorem Kintchine and Levy [ |L7| , 
we express the stable distribution by means of character- 
istic functions with the form exp(i"fk — c\k\ a ). The range 
of a is: < a < 2 , namely wider than that compati- 
ble with the dynamic treatment of this letter that only 
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focuses on the interval: 1 < a < 2. Our treatment is 
restricted to the case of symmetric distribution, thereby 
setting 7 = 0. The same theorem by Kintchine and Levy 
establishes that []: 

c = —aM (a) (ci + c?) cos (^°^ c > 0, 1 < a < 2 , 

(A3) 



where 
M(a) 



(e-y -1+y) 



dy _ T(2 - a) 



yi+'x a(a - 1) 



(A4) 



and Ci and ci are two constants establishing the distri- 
bution asymptotic properties. 

Let us consider the distribution function F(x) of the 
random variable £. According to a theorem established 
by Gncdenko |l7| , the distribution F(x ) co nverges to a 
stable function V(x) in the sense of Eq. ( A2) if, and only 
if, the following relations 



F{x) = {c l a'* + q l {x))^- ifx<0 



and 



F{x) = 1 - (c 2 c 



■Qa(x))— ifx>0 



(A5) 



(A6) 



hold. Note that the function qi(x) and qi (x) vanish for 
x — > — oo and x — > oo , respectively. The constant a 
depends on the choice of B n appearing in Eq.(Al). We 
make the choice B n = n x / a yielding a = 1. Thus in the 
the asymptotic limit \x\ — > oo the distribution F{x) only 
depends on c\ and ci- Furthermore, the earlier choice of 
a symmetric distribution yields c\ — C2- 

Note that to realize our purposes we have to adopt the 
non-normalized form: 



Avi + Av2 



Av n 



(A7) 



It is evident that this non-normalized f orm can be re- 
lated to the normalized form of Eq. (Al) by setting 
Av k = (k n ^ 1 ^ a an d A n = 0. This is equivalent to replac- 
ing the parameter c\, responsible, as earlier shown, for 
the asymptotic distribution properties, with cyn. Claim- 
ing no rigour, we set n = t/T, where T is the mean 
value of the interval between two realizations of Au/c, so 
n is the mean number of realizations of the random vari- 
able Av. Thus we obtain for the stochastic Levy process 
cxp(-6t|fc| Q ): 



By identifying F(x) with the function distribution of 
(|l2|), we obtain: 



(/3TY J+1 



2 (J3T + x/W)P +1 



if x > 0. 



(A9) 



Finall y, b y applying Eq. ( |A6| ) to Eq. (A£) and using 
Eqs. ( A3j) and flAgj) , we obtain the result of Eq. Q). 



[1 

[2 

[3; 

[4 
[5 

[6 

[7; 
[8 

[9 

[10 

[11 
[12 

[13 

[14 

[15 

[16 
[17 

[is; 

[19 
[20 



M.F. Shlesinger, G.M. Zaslavsky and J. Klafter, Nature 
363 (1993) 31. 

J. Klafter, M.F. Shlesinger and G. Zumofen Phys. Today 
49 (1996) 33. 

P. Levy, in: Theorie de l'Addition des Variables Aletories 
(Gauthier-Villars, Paris, 1937). 
J.H. Doob, Ann.Math. 43 (1942) 351. 
E.W. Montroll and B.J. West, in: Fluctuation Phe- 
nomena, E.W. Montroll and J.L. Lebowitz, eds. (North- 
Holland, Amsterdam, 1979). 

V. Seshadri and B.J. West, Proc. Nat. Acad. Sci. USA 
79 (1982) 4501. 

B.J. West and V.Seshadri, Physica A 113 (1982) 203. 
S. Jespersen, R. Metzler, H.C. Fogedby, Phys. Rev. E, 
59 (1999) 2736. 

G.M. Zaslavsky, in: Physics of Chaos in Hamiltonian Sys- 
tem (Imperial College Press, 1998). 
J. Klafter, G. Zumofen, Physica A 196 A (1993) 102. 
G. Zumofen, J. Klafter, Phys. Rev. E 47 (1993) 851. 
G. Trefan, E. Floriani, B.J. West, and P. Grigolini, Phys. 
Rev. E 50 (1994) 2564. 

T. Geisel and S. Thomae, Phys. Rev. Lett. 52 (1984) 
1936. 

P. Allegrini, P. Grigolini and B.J. West, Phys. Rev. E 54 
(1996) 4760. 

T. Geisel, A. Zacherl, and G. Radons, Z. Phys. B 71 
(1988) 117. 

J.Klafter, G. Zumofen, Europhys. Lett. 25 (1994) 565. 
B.V. Gnedenko and A.N. Kolmogorov, in: Limit dis- 
tributions for sums of independent random variables 
(Addison- Wesley, 1954) . 

T. Geisel, J. Heldstab and H. Thomas Z. Phys. B 55 
(1984) 165. 

M. Buiatti, P. Grigolini, L. Palatella, Physica A 268 
(1999) 214. 

M. Bolog na, P. Grigolini, J. Riccardi, submitted to Phys. 
Rev. E, (|cond-mat 9907464|). 



b = c/T 



(A8) 



^ue to a misprint in Eq.(ll) of chap. 7 §34 of Ref. |1J, the 
factor a of Eq. (A3) is missing. 
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